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ABSTRACT 

We study the evolution of the red galaxy fraction (fred) in 905 galaxy groups with 0.15 < z < 0.52. 
The galaxy groups are identified by the 'probability Friends-of-Friends' algorithm from the first Red- 
Sequence Cluster Survey (RCSl) photometric-redshift sample. There is a high degree of uniformity 
in the properties of the red-sequence of the group galaxies, indicating that the luminous red-sequence 
galaxies in the groups are already in place by z 0.5 and that they have a formation epoch of zk,2. 
In general, groups at lower redshifts exhibit larger fred than those at higher redshifts, showing a 
group Butcher-Oemler effect. We investigate the evolution of fred by examining its dependence on 
four parameters, which can be classified as one intrinsic and three environmental: galaxy stellar mass 
(M*), total group stellar mass (M^^g^p, a proxy for group halo mass), normalized group-centric radius 
(Vgrp), and local galaxy density (S5). We find that M, is the dominant parameter such that there is 
a strong correlation between fred and galaxy stellar mass. Furthermore, the dependence of fred on 
the environmental parameters is also a strong function of M». Massive galaxies (M* ^ilO^^ M©) show 



little dependence of fred on r. 
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, and E5 over the redshift range. The dependence of fred on 



these parameters is primarily seen for galaxies with lower masses, especially for M^ <10 Mq. We 
observe an apparent 'group down-sizing' effect, in that galaxies in lower-mass halos, after controlling 
for galaxy stellar mass, have lower fred- We find a dependence of fred on both rgrp and E5 after the 
other parameters are controlled. At a fixed rgrp, there is a significant dependence of fred on S5, while 
fgrp gradients of fred are seen for galaxies in similar S5 regions. This indicates that galaxy group 
environment has a residual effect over that of local galaxy density (or vice versa), and both parameters 
need to be considered. This result suggests that processes identified with local galaxy density, such 
as galaxy harassment and mergers, and those associated with accretion into a larger group halo, such 
as ram pressure and strangulation, are both partaking in driving galaxies to their final red quiescent 
state. We discuss these results in the context of the 'nature vs nurture' scenario of galaxy evolution. 
Subject headings: galaxies: evolution - galaxies: photometry - galaxies; cluster: general 



1. INTRODUCTION 

Galaxies evolve with time. Generally speaking, the 
star formation rate within galaxies decreases as galaxies 
age, and hence galaxy colors transit from blue to red, 
and their spiral arms become less and less dominant. It 
has been known since the 1970's that red galaxies tend to 
populate galaxy clusters and blue galaxies are common 
in the field (Oemler 1974). The fraction of red galax- 
ies decreases from cluster center to the field. Measuring 
galaxy morphologies in 55 clusters in the nearby Uni- 
verse, Dressier (1980) found that galaxy population frac- 
tions have a stronger dependence on local galaxy density 
than on cluster-centric radius. This work formulated the 
well-known 'morphology-density' relation, and is gener- 
ally investigated as the 'color-density' relation later on 
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in the literature (e.g., Kodama et al. 2001; Cooper et al. 
2006; Quadri et al. 2007; Cassata et al. 2007; Cooper et 
al. 2007). This result is interpreted as a decrease in the 
star formation rate as local galaxy density increases. A 
'critical density' was suggested to characterize the local 
galaxy density effect (e.g., Lewis et al. 2002; Gomez et 
al. 2003; Tanaka et al. 2004) with the concept that star 
formation rate and galaxy color exhibit abrupt changes 
when the galaxy density crosses the critical density. Re- 
gardless of whether this critical density exists or what its 
true meaning is, the density where both star formation 
rates and galaxy colors change rapidly may indicate an 
environmental transition between galaxy groups and the 
field (e.g., Kodama et al. 2001; Goto et al. 2003; Balogh 
et al. 2004; Li et al. 2009, hereafter L09). 

Galaxy groups are small galaxy aggregations bound by 
gravity. Some groups contain a few galaxies while oth- 
ers may have several dozens. This gives galaxy groups 
a fundamental role in building up large-scale structures, 
as more and more galaxies are accreted into groups over 
time and groups merge together. In the local Universe 
the majority of galaxies are found in galaxy groups (e.g., 
Geller & Huchra 1983; Fukugita et al. 1998; Eke et al. 
2004). Because of their sufficiently high galaxy density 
and low velocity dispersion, galaxy groups are a favored 
environment for interactions and mergers, causing most 
of the galaxy transformations in morphology, star forma- 
tion rates, and colors (e.g., Ghigna et al. 1998; Mulchaey 
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& Zabludoff 1998; Hashimoto & Ocmler 2000). The sup- 
pressed star formation in galaxy groups hence will lead to 
a decline in the mean star formation rate of the Universe 
as structures grow with decreasing redshift. This sce- 
nario provides a plausible explanation for the declining 
global star formation rate with cosmic time (e.g., Madau 
et al. 1996; Lilly et al. 1996; Cowie et al. 1999; Cooper 
et al. 2008). 

Within a galaxy group, the properties of galaxies ex- 
hibit a dependence on their distance from the group cen- 
ter. Observationally, galaxies in the center of groups 
are more likely to be bright, red, and early-type. The 
fraction of early-type galaxies decreases with increasing 
group-centric radius, while the fraction of fainter, bluer, 
and late-type galaxies has a positive correlation with 
group-centric radius (e.g., Dommguez et al. 2002; Brough 
et al. 2006; Weinmann et al. 2006). Galaxies in the cen- 
ters of low-mass groups already show different properties 
from those in the field (e.g., Balogh et al. 2007; Kawata & 
Mulchaey 2008). Such differences are stronger in higher- 
mass groups, and galaxies in these groups have prop- 
erties close to those in clusters (e.g., Domfnguez et al. 
2002). Galaxy groups, therefore, can be considered as 
'mini-clusters'. Yet galaxy groups are systems typically 
smaller than clusters by a factor of ^-^10 or more in mass. 
Mechanisms such as ram-pressure stripping and harass- 
ment are perhaps not important within galaxy groups. 
This makes galaxy groups useful sites to study galaxy 
evolution possibly free of complicated multiple driving 
mechanisms. 

Most studies on galaxy groups, however, are in the lo- 
cal Universe using samples from the Sloan Digital Sky 
Survey (SDSS) and the Two-degree-field (2dF) Survey 
(e.g., Balogh et al. 2004; Martinez & Muriel 2006; Pog- 
gianti et al. 2006; Robotham et al. 2006). At ^ > 0.1, our 
current understanding of galaxy groups is based on sam- 
ples considerably smaller than that available in the local 
Universe; for examples, the CN0C2 (e.g., Carlberg et al. 
2001; Wilman et al. 2005; Wilman et al. 2005; Balogh et 
al. 2009), DEEP2 (e.g.. Cooper et al. 2007; Gerke et al. 
2007), and zCOSMOS (e.g., Knobel et al. 2009; lovino 
et al. 2010) group catalogs. Galaxy groups are not easily 
detected beyond the local Universe because of their low 
gas content, as well as the less significant number-density 
contrast to the background. Traditionally, galaxy groups 
are identified within spectroscopic-redshift samples using 
algorithms such as the friends-of-friends or the Voronoi 
tessellation method (e.g.. Carlberg et al. 2001; Gerke et 
al. 2005; Berlind et al. 2006; Sochting et al. 2006; Tago 
et al. 2008). The requirement of large spectroscopic sur- 
veys with high completeness adds difficulties to the study 
of galaxy groups at ^; > 0.1, since obtaining such cata- 
logs requires huge amounts of time, money, and work. 
Another route to identify galaxy groups is to apply ded- 
icated group-finding algorithms to photometric-rcdshift 
samples (e.g., Botzler et al. 2004; Li & Yee 2008). Galaxy 
groups found this way have much larger uncertainties in 
their redshift estimation, and contain background galax- 
ies as their members. The contamination of false groups 
may also bias scientific analyses. Even so, the ability 
to obtain large samples of photometric-rcdshift galaxy 
groups will still allow us to gain useful information, after 
correcting for background contamination statistically. 

We study galaxy groups at 0.15 < z < 0.52 using a 



photometric-rcdshift group sample drawn from the multi- 
band photometric data of the first Red-Sequence Cluster 
Survey (RCSl; Gladders & Yee 2005). We aim to explore 
how colors of galaxies, which are a proxy for galaxy popu- 
lation, change with redshift, and how environment affects 
galaxy populations therein. We follow the methodology 
used in our study of groups associated with CNOCl clus- 
ters (L09) in deriving photometric redshifts, identifying 
galaxy groups, and defining environmental parameters. 
The structure of this paper is as follows. We describe 
our data and group sample in §2 and §3. The results 
are presented in §4, and discussed in §5. We summarize 
our work in §6. We adopt the cosmological parameters 
i?o=70 km/s/Mpc, O„=0.3, and Oa=0.7. 

2. THE DATA 
2.1. The Observation 

We derive our galaxy group sample from the Northern 
patches of the RCSl, which is an imaging survey with 
the primary goal of measuring fl„i and (Ts using the clus- 
ter mass function (Gladders & Yee 2005; Gladders et al. 
2007). The survey covers 22 widely separated patches 
in total in the Northern and Southern sky with a total 
area of 92 square degrees. The observations were carried 
out from May 1999 to January 2001. The imaging in i?c 
and z' of the ten Northern patches was conducted using 
the CFH-12K camera on the 3.6m CFHT telescope. The 
camera has 12kx8k pixels in total, consisting of twelve 
2kx4k CCDs with 0.206 " per pixel, providing a 42' x 28' 
field of view. The layout of each patch is arranged with 
15 pointings in a slightly overlapping grid of 3 x 5 point- 
ings, giving a patch size of 2.1 x 2.3 deg^. The integra- 
tion times are 900 and 1200 seconds with average seeing 
of ^ 0.70" and ^ 0.62" in the Rc and z' passbands, re- 
spectively. The 5cr limiting magnitudes for point sources 
are Rc = 24.8 (Vega) and zab = 23.9 in an aperture of 
diameter 2.7". The details of the data reduction and pho- 
tometric catalog can be found in Gladders & Yee (2005). 

To obtain photometric redshifts, we require additional 
photometric bands. The B and V photometry for the 
RCSl Northern patches was obtained as a follow-up 
project using the CFH-12K camera (Hsieh et al. 2005). 
It covers 33.6 deg^ in total with 108 pointings (^ 75% of 
the original RCSl-CFHT patches). This defines the sky 
area used for our group catalog. The runs were carried- 
out from May 2001 to June 2002. The typical exposure 
times in B and V are 840s and 480s, respectively. The 
average seeing is ~ 0.95" in B and ~ 0.65" in V. The 
5(7 limiting magnitudes are B=25.0 and V=24.5, on av- 
erage, within an aperture of 2.7" diameter. We refer to 
Hsieh et al. (2005) for more details on the follow-up data 
set. 

2.2. The Photometric- Redshift Catalog 

We use the empirical photometric-redshift method 
modified from Li & Yee (2008) to estimate galaxy red- 
shifts. This photometric redshift algorithm assumes that 
galaxy redshift is a polynomial function of galaxy mag- 
nitudes and colors, and the coefficients are derived from 
a training set which is a catalog containing spectroscopic 
galaxy redshifts and multi-band photometry. We con- 
struct the redshift training set using spectroscopic sam- 
ples from the Hubble Deep Field (HDF; Giavalisco et 
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al. 2004), the Canadian Network for Observational Cos- 
mology Survey (CN0C2; Yee et al. 2000), and the Deep 
Extragalactic Exploratory Probe (DEEP2; Vogt et al. 
2005). The photometry of the training set is processed 
and calibrated in the same manner as the RCSl catalogs, 
and the details can be found in Hsieh et al. (2005). The 
final training set contains ^5,300 galaxies with BVRz' 
photometry. 

Instead of dividing the training set and all input galax- 
ies into several fixed color-magnitude cells as was done 
in Li & Yee (2008), we derive the coefficients of the 
photometric-redshift polynomial fit individually for each 
input galaxy by using a subset of the training set which 
contains 400 galaxies chosen from the complete training 
set whose magnitudes and colors are the closest to the 
input galaxy. The choice of the size of the sub-training 
set is a compromise between accuracy and computation 
time. Our tests show that any sub-training set with 
more than 250 galaxies will provide reasonable photo- 
metric redshift results. These training galaxies are cho- 
sen based on quadratically summed ranks of color and 
magnitude differences between the training set galaxies 
and the input galaxy. All four magnitudes and six colors 
are used. In doing the fitting, we further assign weights 
to the chosen training-set galaxies, based on the inverse 
value of their final rank. This partially alleviates the 
redshift bias that may be introduced to input galaxies 
near the edges of the color-magnitude distributions. The 
photometric-redshift errors are estimated empirically by 
assuming Gaussian magnitude errors and bootstrapping 
the training set in the color-magnitude cell. We remove 
galaxies with very large photometric redshift uncertain- 
ties from the sample. A galaxy is considered not to have 
an acceptable photometric redshift if its uncertainty is 
larger than the very loose criterion of 0.6(1 +z). As in 
L09, we assign a weight, Wi, to each galaxy based on 
the inverse of the fraction of galaxies with acceptable 
photometric redshifts to the total as a function of Rc 
magnitude. This Wi is found not strongly dependent on 
galaxies colors (Yee et al. 2005; Li et al. 2009). The effect 
of Wi on our results will be further discussed in §5.1. 

To test the photometric-redshift accuracy, we ran- 
domly exclude 200 galaxies from the training set and 
treat them as input galaxies. Their photometric red- 
shifts are then derived using the rest of the training set 
galaxies, and compared with their known spectroscopic 
redshifts, as presented in Figure 1, along with their es- 
timated uncertainties. For galaxies within our redshift 
range of interest (0.15 < Zgpec < 0.6), the overall la dis- 
persion (Figure 1) is ^-^0.064 with a mean Zp^ot — Zgpec 
off set of 0.011. More tests of our photometric-redshift 
technique using both simulated data and the GOODS-N 
catalog can be found in Hildebrandt et al. (2010). 



3. 
3.1. 



THE GALAXY GROUP SAMPLE 



Group Catalogs and the Sample 

The identification of galaxy groups is achieved using 
the 'probability Friends- of- Friends' (pFoF) algorithm, 
described and tested in detail in Li & Yee (2008). The 
algorithm uses the group redshift probability density as 
a conditional probability density and the photometric- 
redshift probability distribution of individual galaxies to 
perform the friends-of-friends linkage in the third dimen- 
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Fig. 1. — Panels a) and c): Comparisons between photometric 
and spectroscopic redshifts for 200 galaxies whose photometric red- 
shifts are derived using the rest of the ~5100 training set galaxies. 
The \a dispersion of Zpf^^t — Zspec is ~0.064 for 0.15 < z < 0.6. 
Panel b): Empirically estimated l-cr photometric redshift uncer- 
tainties of individual galaxies as a function of Zspec- 

sion. A major feature of this algorithm is that the group 
redshift probability density is improved as more linked 
group members are added. 

The pFoF algorithm has been tested using a simulated 
catalog (Croton et al. 2006) based on the Virgo Con- 
sortium Millennium Simulation (Springel et al. 2005). 
Different linking criteria in both the 2D sky plane and 
the redshift direction, as well as various sample depths, 
have been adopted to investigate the group-finding per- 
formance. The results of these tests are described at 
length in Li & Yee (2008). The RCS group catalogs are 
obtained by applying the pFoF algorithm to a sample of 
depth M^^ -f 2.0, where M^^ =-21.41 (Kodama & Ari- 
moto 1997). We have generated a complete catalog of 
galaxy groups which is available in electronic form (Ta- 
ble 1). Galaxy groups in this catalog are selected to have 
Ngai > 5 and Ng^ > 5, and redshift between z—Q.15 and 
z=Zcut- Here, Ngai is the net weighted member count to 
Af^^+2.0 (i.e., Ngai = J2i''^i) ~ -^bg , where the computa- 
tion of the background counts N^g is presented in §3.5); 
Ngz is the actual number count of the linked galaxies in 
the pFoF group; and Zcut is the redshift where the nomi- 
nal Rc limit of a patch has a depth of M'^ -1-1.5. The Rc 
limit is defined as the magnitude at which Wi reaches 2.0, 
and ranges from 21.16 to 22.86 for the different patches, 
corresponding to a Zcut range of 0.30 to 0.53. The Ngai 
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and Ngz limits of > 5 arc chosen as a compromise be- 
tween having a sample reaching down to sufBciently low 
group mass and the number of false detections, based on 
the results of our mock tests. 

These criteria produce a sample of 1153 groups. Be- 
cause we use a photometric redshift technique modified 
from that of Li & Yee (2008), wc redo the mock tests 
for the pFoF algorithm. The results, which are similar 
to those in Li & Yee (2008), show that this sample may 
contain ~35% false groups (see §5.1). However, we can 
use other information to attempt to cut down the frac- 
tion of false groups. For examples, the total stellar mass 
of a group INL grp and the richness parameter Bgc are 
measured using net and galaxy counts, respectively, 
within an aperture on the sky, rather than based on FoF 
connected galaxy counts, and thus provide different infor- 
mation (see §3.4 for the calculations of M* and Bgc). 
Thus, groups with small M^.g^p or Bgc, and in particular 
groups with very discrepant M^ g^pand Bg,. (i.e., large 
BgcSvaaXl M^^g^p, or vice versa) are expected to more 
likely be false groups. This is verified by our mock group 
sample tests. Applying the same M*,grp cut to the mea- 
sured M, grp in our mock catalog, we find that ~ 65% 
of the false groups have log(M*,grp/M0)<11.2. For our 
analyses, wc further employ a total group stellar mass cut 
of log(M,_grp/M0)>11.2 and a richness cut of Bgc > 125 
Mpc^-^/igQ^'* in selecting our group sample. The total 
group stellar mass and richness cuts remove 248 groups 
from the sample. Assuming that these arc mostly false 
detections, we expect ^-^15% false groups in the remaining 
group catalog, which we use as our final group sample. 
We note that the most significant effect that false groups 
produce in our analysis is to bias the galaxy red fraction 
towards smaller values (sec §5.1 and Li & Yee (2008)), 
especially for samples of poor groiips where the false de- 
tection rate is higher. 

We present the basic data of 24 galaxy groups in Table 
1 as examples. They are typical groups in each redshift 
and M^^grp bins in our analyses in §4. Figures 2, 3, and 
4 show the sky locations and observed color-magnitude 
diagrams of these groups. The full sample of 1153 groups 
is presented in the electronic version of the Table. Some 
groups are also found as RCS clusters (Gladders & Yee 
2005). We note that most groups, especially the richer 
ones, exhibit a clear red sequence at the expected theo- 
retical model colors. 

The redshift distribution of our group sample is plotted 
in Figure 5. Since different patches have different redshift 
limits, we also plot in Figure 5 the redshift distribution 
corrected for the sampling areas for the different redshift 
bins. For our analysis, we divide the groups into three 
redshift bins: 0.15 < z < 0.35, 0.35 < z < 0.45, and 
0.45 < z < 0.52 designated as the z 0.25, z ^ 0.4, and 
z ~ 0.5 bins in our analysis. There are 304, 317, and 284 
groups, respectively, in these bins. 

3.2. Samples of Group Galaxies 

To select galaxies in the same redshift space as a group, 
we follow the method used by L09 for selecting galaxies 
in the same redshift space as the CNOCl clusters. That 
is, the photometric-redshift probability for a galaxy to 
be at the same redshift as a given group has to satisfy 
a threshold, which is set to be the same as the redshift 
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Fig. 2. — Examples of sky locations and observed B — Rc versus 
Rc color-magnitude diagrams for RCS galaxy groups. The group 
ID and redshift arc indicated in the title of each plot. The circles in 
the sky maps have a radius of one -R200 • The linked pFoF members 
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TABLE 1 
Catalogs of RCS Galaxy Groups" 



ID 



R.A.' 



Dec' 



M' 



1447_48 

2148.59 

1120.57 

0920.317 

1417.181 

1614.469 

2316.15 

2148.153 

0223.133 

1417.23 

0920.55 

0920.121 

0920.365 

1614.70 

0920.174 

0920.322 

1614.54 

2316.41 

1614.562 

1417.321 

1447.45 

0920.28 

0223.252 

0223.246 



14:47:59.4 
21:52:33.6 
11:22:17.1 
09:22:21.9 
14:16:54.0 
16:14:36.1 
23:14:15.5 
21:53:54.8 
02:27:53.7 
14:14:35.6 
09:23:0.31 
09:31:23.5 
09:29:2.96 
16:13:30.8 
09:22:0.56 
09:28:5.15 
16:17:18.1 
23:16:59.4 
16:12:41.3 
14:19:0.69 
14:46:45.5 
09:23:42.5 
02:22:19.8 
02:25:57.5 



09:52:46.0 
-05:29:24.9 
25:58:6.66 
37:50:15.7 
52:21:29.4 
30:25:51.2 
00:43:42.8 
-05:15:33.7 
01:14:10.1 
53:54:13.4 
37:32:48.4 
37:50:9.19 
37:49:43.8 
30:00:6.82 
37:54:28.6 
37:55:16.3 
30:28:27.9 
-00:12:32.5 
29:55:42.9 
53:51:41.1 
10:07:3.18 
37:20:13.6 
00:05:12.9 
00:54:40.9 



0.198±0.011 
0.277±0.015 
0.349±0.019 
0.375±0.010 
0.487±0.018 
0.490±0.014 
0.405±0.012 
0.443±0.015 
0.207±0.010 
0.253±0.008 
0.380±0.010 
0.503±0.018 
0.497±0.017 
0.365±0.012 
0.418±0.011 
0.430±0.017 
0.400±0.009 
0.442±0.016 
0.491±0.019 
0.519±0.023 
0.263±0.010 
0.220±0.009 
0.390±0.011 
0.396±0.016 



6.014 
7.947 
8.082 
5.486 
6.117 
6.512 
11.18 
6.341 
5.961 
6.497 
8.591 
9.267 
7.045 
8.067 
6.094 
5.371 
7.978 
10.36 
5.272 
6.664 
6.237 
5.666 
5.754 
5.206 



12.45 
12.45 
12.41 
12.41 
12.41 
12.41 
12.39 
12.39 
12.02 
12.02 
12.01 
12.01 
12.01 
12.01 
11.96 
11.96 
11.65 
11.64 
11.62 
11.61 
11.58 
11.57 
11.56 
11.56 



640±211 
271±174 
442±211 
653±229 
674±250 
167±203 
685±233 
454it223 
736±218 
638±205 
808±243 
244±218 
213±213 
443±205 
344±208 
510±226 
366±203 
395±214 
207±207 
259±223 
188±173 
356±175 
334±205 
264±200 



''The photometric catalogs are available in Hsieh et al. (2005) 
''in J2000. 

'^Ngai- net pFoF member count after correcting for background 
contamination and completeness weights. 

Mgrp=log{Mii:.grp/MQ), computed using all group members 
within 0.5-R200 and with background contamination corrected. 

^Bgc: computed using all galaxies within 0.25Mpc to the group 
centers 

linking criterion in the pFoF algorithm. These galaxies 
are called 'group galaxies' in our analyses (§4) and are 
not limited only to the galaxies linked by the pFoF algo- 
rithm, since redshift is the only criterion used for their 
inclusion. In essence, the main task of of the pFoF algo- 
rithm is to identify the locations where there arc groups; 
the use of the linkage criterion in redshift then gener- 
ates a sample of galaxies whose photometric redshifts 
are consistent with being in these groups. In our anal- 
yses, we further impose group-centric radius criteria for 
the selection of group galaxies in various radial bins. We 
have in total 5028, 5832, and 5560 group galaxies (before 
background subtraction) brighter than -f 1.5 within 
0.5ii!200 (derived from group richness, see §3.3) in the 
three redshift bins. 

We also estimate the stellar mass for each group galaxy 
using an empirical relation: log(M,/ii^^) — —0.523 + 
0.683(Mb - Mr J, from Bell et al. (2003) with solar 
units from Worthey (1994). The Mb and Mr^ are the k- 
corrected rest-frame absolute magnitudes derived using 
the same method as L09. The k-correction values used 
for each galaxy are based on its color and photometric 
redshift. 



Figure 6 presents the plot of absolute magnitude ver- 
sus stellar mass for our group galaxy sample. The galaxy 
sample is complete to M^^ -|- 1.5 in all three redshift 
bins. In stellar mass, the sample is essentially complete 
to log(M,/M0)=lO.2 for the z ~ 0.25 and z ~ 0.4 
samples; while the z ^ 0.5 sample is 100% complete 
only to log(M*/M0) ~10.8. We note that the level 
of incompleteness at z ^ 0.5 is corrected substantially 
by the Wi weight assigned to each galaxy. There are 
a total of 3550, 3933, and 3940 group galaxies within 
0.5i?2oo with log(M*/M0)>lO.2 in the 3 redshift bins, 
respectively. The number counts for group galaxies at 
0.5 < i?2oo < 10 are 2637, 3298, and 4605 galaxies for 
the same redshift divisions. 

3.3. Group-Centric Radius and Scaling Radius 

To define the center of a pFoF group, we first take the 
median R.A. and Dec. of all the linked members as the 

'fiducial' group center. Each linked galaxy is assigned a 
score based on its Rc magnitude and distance to the fidu- 
cial center. The score is then computed as the quadratic 
sum of the ranks of the magnitude and distance differ- 
ences between the fiducial group center and the galaxy. 
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Fig. 3. — Same as Figure 2. 



Fig. 4. — Same as Figure 2. 
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Fig. 5. — The redshift distribution of 905 galaxy groups in our 
sample using a 0.01 bin size in z. The thin red and thick black 
histograms are the distributions before and after applying area 
corrections. 
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Fig. 6. — Stellar mass M* and M^'^ for group galaxies within 
0.5iJ200 from 182 groups in the RCS 0920 patch. The ver- 
tical dashed line marks Mt. + 1.5 while the horizontal one, 
log(M*/M0)=lO.2. The red and blue solid lines mark the ridge 
lines for red and blue galaxies, respectively. Selecting the sample 
using stellar mass rejects many blue galaxies. 

Both the magnitude and distance contribute the same 
weight to the score. The galaxy which has the best score 
is defined as the 'central galaxy'. We adopt the position 
of the central galaxy as the final center position of the 
group from which a group-centric radius is defined. If 
more than one linked galaxy in the same group has the 
same score, we choose the center as the one with the 
brighter magnitude. 

For each group galaxy, we compute the projected dis- 
tance to the center and scale it by R20Q, the radius within 
which the density is 200 times the critical density, so that 
our analysis is not limited only to either the cores of rich 
groups or the outskirts of poor ones. The i?2oo of each 
group is estimated from the correlation between cluster 
richness Bgc (see §3.4) and i?200 derived using the X- 
ray luminous CNOCl clusters in Yee & Ellingson (2003) 
(also see Barkhouse et al. (2007)). The typical i?200 in 
our RCS group sample is ^ 1.07 ± 0.24Mpc. The group- 
centric radius is denoted by r^^p in units of i?200- We 
note that the Bgc — i?200 relation for low-mass galaxy 
groups may not be identical to that of the CNOCl clus- 
ters, but it still gives us an estimate and rank-order of 
i?2oo- The effects of the uncertainties in the i?200 on our 
results are further discussed in §5.1.3. 

3.4. Group Richness and Masses 

Galaxy groups and clusters with more members are 
likely associated with more massive dark matter halos 




12.0 12.5 

l09(M..„p/Me) 

Fig. 7. — The histogram of total stellar mass within 0.5iJ200 in 
each galaxy group. The vertical dashed line is the median value, 
and the two (red) dotted lines are the 25 and 75 percentiles. The 
median total group stellar mass remains similar over the three red- 
shift bins. 

and have deeper gravitational potentials (e.g., Yee & 
Ellingson 2003; Popesso et al. 2007). Therefore, to inves- 
tigate the influence caused by galaxy groups, the group 
richness should be considered. We explore the param- 
eterization of group richness using two different param- 
eters: the total stellar mass within a fixed fraction of 
i?200 and the richness parameter Bgc (Longair & Seldner 
1979). 

We compute the total stellar mass of a galaxy 
group, M*,grp, using the sum of the stellar masses 
of all the group galaxies within 0.5 i?200 in the log 
(M*/Afs„„)>10.2 sample. The M^,^grp is corrected for 
background contamination (see section §3.5). We use 
0.5i?2oo instead of li?2oo because group member counts 
have a larger excess relative to the background within 
the smaller radius. We note that group galaxies in the 
z ~ 0.5 bin are incomplete at log(M*/AfQ)=10.2. How- 
ever, partly due to the galaxy weight corrections and the 
small contribution to the total group mass from low-mass 
galaxies, this has only a few percent effect on the total 
group stellar mass. This is found to be the case when the 
incompleteness effect is tested using the average ratios of 
M*_grp computed with limits of log(M*/M0)=lO.2, 10.6, 
and 11.0 for each redshift bin. Figure 7 presents the 
^*,grp distribution of our group sample in each redshift 
bin. The distributions are similar for the three redshift 
bins, peaking at or slightly above log(M*/M0)=12.O, 
with median values differing by ^0.1 dex. The distri- 
butions also show that the sample is beginning to be 
incomplete at log{M^, ,grp / Mq) s12. This is likely due to 
the pFoF algorithm's inability to find loose poor groups, 
and conclusions based on analyses of the low-mass groups 
should take this bias into consideration. 

The Bgc parameter is defined as the amplitude of the 
cluster center-galaxy correlation function, and can be 
estimated via a deprojection of the angular correlation 
function onto the spatial one assuming spherical symme- 
try (Longair & Seldner 1979). It is estimated by counting 
excess galaxies up to a certain absolute magnitude within 
a given radius, and is corrected for the background counts 
and scaled by an average luminosity function and average 
spatial profile. Details on its computation, robustness as 
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a cluster richness parameter, and as a mass indicator are 
presented in Yee & Lopez-Cruz (1999) and Yee & EUing- 
son (2003). 

For the RCS-1 clusters, Bgc is computed in a 0.50 Mpc 
h^Q radius (Gladders & Yee 2005). Considering that 
galaxy groups are poorer systems than galaxy clusters, 
we calculate Bgc for each galaxy groups using a 0.25 Mpc 
h'^Q radius to minimize background noise. We find that 
the Bgc values computed using a radius of 0.50Mpc are 
~ 35% smaller than those using a 0.25Mpc radius, in- 
dicating a steeper slope of the galaxy correlation func- 
tion than the canonical 7 = 1.8. We note that a cluster 
of Bgc = 600Mpci-S/i^o^-^ (within a radius of 0.5Mpc) 
has a richness equivalent to Abell class 0, or a mass of 
~ 3 X 10"/i-iMo (Yee & Ellingson 2003). 

Figure 8 presents the comparison between Bgc and 
^*,grp- A positive correlation is observed, with the scat- 
ter broadly consistent with the measurement errors. We 
obtain the best fitting power-law, using the more robust 
Ngai > 8 sample: log(M,,g^p/MQ) lAl\og{Bgc) + 8.30. 
Using the simple scaling relations in Yee & Ellingson 
(2003) with the assumption that galaxies trace total halo 
mass, we should expect the halo mass within a dynam- 
ically scaled radius (e.g., M20Q, the mass within i?20o) 

3/'Y 

to scale as Bgc , where 7 is the power-law slope of the 
galaxy-galaxy correlation function; or, for the canoni- 
cal 7 of 1.8, log(M2oo) '~ Bgc^- Thus, our power-law 
fit of M*_grp and Bgc is consistent with M^^grp being 
a direct one-to-one proxy for the Af2oo halo mass of 
the groups. Using the correlation from Yee & Elling- 
son (2003) between Bgc and M200, and compensating 
for Hq, we can roughly estimate the relation between 
the total group stellar mass and the halo mass. We do 
this by fitting the Mst^grp^Bgc relation using the same 
power-law slope (1.64) as that measured by Yee & Elling- 
son (2003) for the M2oo-Bgc relation for the CNOCl 
clusters, and correcting for factor of 0.65 in Bgc val- 
ues computed using a 0.5 /i^g^ Mpc radius. This gives 
log(M2oo,grp)'^log(M*^grp)+l-85. This provides a very 
rough scaling, probably no better than a factor of two, 
between the group stellar mass and M200. In this scaling, 
a group with log(M,^grp/Afo)=12.0, (approximately the 
median M^^^^p for our sample), represents a halo mass of 
~ 7 X lO^'^M©. In the remainder of the paper, we use the 
group stellar mass M+^g^p as the richness indicator and a 
rough halo mass proxy. 

3.5. Background Correction 

Because of the relatively large uncertainty of the pho- 
tometric redshift, even with the use of the pFoF algo- 
rithm for determining whether a galaxy belongs to a 
given group, there will still be a significant fraction of 
"group galaxies" which are field, or near-field, galaxies, 
especially at the larger group-centric radii. Hence, it is 
essential to perform background contamination correc- 
tions in analyses using the group galaxy sample. We con- 
struct both luminosity-limited and stellar-mass-selected 
samples of background counts using three RCS patches: 
0920 , 1417 ,and I6I4 ■ These three patches have the 
deepest photometry and cover a total of 11.38 deg^ on 
the sky. We note that these three control patches are 
part of our whole sample. The deep photometry allows 
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Fig. 8. — The comparison between M«,grp and Bgc, two parame- 
ters used as indicators of group richness. A correlation is seen. The 
dotted lines indicate the Bgc=125 and log(M*,9rp/MQ)=11.2 crite- 
ria for selecting our final group sample. The solid line is the linear 
fit (displayed on the top of the figure) between log(M»_grp/M0) 
and log(iJgc) using groups with Ng^i > 8, while the dashed line is 
the fitted relation using a slope of 1.64 (the slope obtained by Yee & 
Ellingson (2003) fitting Af2oo vs Bgc for a sample of CNOCl clus- 
ters). The vertical dashed line marks the mass of M200 ~ 3 X 10^* 
M0, based on the calibration of Yee & Ellingson. Galaxy groups 
are represented by circles, with larger symbols representing groups 
at lower redshift. Squares mark groups with linked pFoF mem- 
bers Ngai > 8. The error-bars on the lower right represent typical 
uncertainties in the two parameters. 

US to compute galaxy counts to a higher redshift. 

We use the method in L09 to construct the galaxy num- 
ber surface density per Mpc^ in the field in redshift bins 
of Az = 0.01, denoted as n(z, Affl^) andn(z,M*). Briefiy, 
the n(z) functions are obtained by integrating the sum of 
weighted photometric redshift probability densities of all 
galaxies in the three patches within each small redshift 
interval and normalized by the area of the patches. For 
a given set of galaxies, e.g., a sample of group galaxies 
within certain group-centric radius over certain magni- 
tude range, the total number of background galaxies is 
computed by multiplying the appropriate n{z) curve with 
the summed photometric-redshift likelihood function of 
all galaxies of interest, normalized by the area. A more 
detailed description of the method can be found in L09. 

3.6. Local Galaxy Density 

It has been shown that local galaxy density plays a role 
in transforming star-forming galaxies into passive ones 
(e.g., Gomez et al. 2003; Goto et al. 2003). We compute 
local galaxy density S5 based on the nearest 5"* neighbor 
galaxy for each group galaxy using the 'group galaxies' 
in our two samples. The S5 is then corrected for back- 
ground galaxy contamination (see L09 for more details 
of the computation). The E5 is computed using galax- 
ies brighter than Af]^ -I- 1.5 for the luminosity-limited 
sample, and using galaxies with log(Af*/M0) > 10.2 
for the stellar-mass-selected sample. We note that in 
general blue galaxies in the stellar-mass-selected sample 
have somewhat smaller E5 than those in the luminosity- 
limited sample, and the differences are primarily seen 
in low E5 regions. This is because most blue galaxies 
in the luminosity-limited sample have stellar mass less 
than log(Af*/M0) = 10.2, and they are excluded in the 
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stellar-mass-selected sample. 

3.7. Red Galaxy Fraction 

As the red sequence of early-type galaxies is the fi- 
nal destination for galaxies on a color-magnitude dia- 
gram, we probe the dependence of the red galaxy frac- 
tion, /red, with various group properties. The fred is 
computed based on the same method as L09 but using 
samples of galaxies defined by stellar-mass limits. The 
red galaxies are defined as those whose colors are red- 
der than half the color differences between the spectral 
energy distributions (SEDs) of E/SO and Sc galaxies at 
a given redshift, which is tabulated in Poggianti (1997). 
This color division changes with redshift on an observed 
color-magnitude diagram, and follows the passive evo- 
lution of galaxy colors. The color division is typically 
0.25 to 0.30 mag to the blue of the red sequence, and 
is comparable to the AB — V = 0.2 used by Butcher 
& Oemler (1984), and the definition used by Ellingson 
et al. (2001). The background count corrections for the 
red galaxy counts are derived using the method described 
in §3.5 but applied to background galaxies with the color 
cut described above. We compute fred using a statistical 
inference (D'Agostini 2004; Andreon et al. 2006) because 
the background contamination plays a role in estimating 
the true fraction of red group members especially when 
the number of galaxies in a group is not large. 

4. RESULTS 

In the following subsections, we investigate fred as a 
function of redshift, M^^g^p, ^grp, and S5. We have ex- 
plored results using both the luminosity and stellar-mass 
samples. We find that both samples exhibit consistent 
trends and lead to similar conclusions. We therefore 
present and discuss our results using only the stellar- 
mass selected sample. We note that, as discussed in 
L09, using a luminosity-selected galaxy sample is equiv- 
alent to having different stellar-mass limits for the red 
and blue galaxies. Our luminosity-selected sample has 
equivalent stellar-mass limits of log(M*/M0)~ 10.5 and 
log(M,/M0)~ 9.8, for red and blue galaxies (Figure 6), 
respectively. On the other hand, the nominal stellar- 
mass limit of log(MH./M0)=lO.2 of the stellar-mass se- 
lected sample corresponds to Mj^'^ ~ —19.0 and -20.75, 
for red and blue galaxies, respectively. Thus, for the 
z ^ 0.5 bin the red group galaxy sample is substantially 
incomplete at log(M»/Af0)=lO.2, and we include results 
from the 2; ~ 0.5 bin only for samples with a stellar- mass 
limit of log(M,/Mo)>10.6. 

To investigate the dependence of galaxy population 
on group richness or mass, we also divide our group 
sample into high-, intermediate- and low-mass bins of 
log(M,,grp/Af0) > 12.2, 12.2>log(M,,g,.p/M0)>11.8, 
and 11.8>log(M,_grp/A^0)>ll-2. These IdIus correspond 
to, approximately, M200 > LIxIO^^^Mq, LIxIO^'^Mq > 
M200 > 4.5 X 10^^ Mq, and 4.5 x lO^^Mg > M200 > 
0.7 x 10^'^ M0, respectively. We intend to probe fred out 
to Vgrp 2. Due to the less accurate redshift informa- 
tion in photometric redshifts, a small number of galaxies 
may belong to more than one group. In this situation, 
we re-assign the galaxy's group membership to the group 
for which the group-centric radius is the smallest. 
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Fig. 9. — Stacked rest-frame B — R vs Mn^ for group galaxies 
within 0.5-R2OO in each redshift and Mt,grp bin. The gray dashed 
hnes in each panel mark log(M, /M0)=11.O, 10.6, and 10.2 from 
left to right. The vertical green line is Mfi^=21Al. Red dots rep- 
resent galaxies satisfying the definition of red galaxies as described 
in §3.7. The solid lines are the fitted red sequences using red galax- 
ies with log(M*/M0)>lO.6. The fitting results are listed in Table 
2. Group galaxies, including those in poor groups, exhibit similar 
color-magnitude relations as those in clusters. 

4.1. Color-magnitude Diagrams and the Red-Sequence 
for Group Galaxies 

Figure 9 presents the stacked rest-frame color- 
magnitude diagrams for group galaxies within 0.5i?2oo at 
each redshift and M^^ g^p bin. We do not apply any back- 
ground correction in Figure 9, but the contamination of 
background galaxies is not expected to be large for rea- 
sonably luminous galaxies, since these group galaxies are 
already pre-selected by photometric redshift and a rela- 
tively small group-centric radius is used. The very small 
number of galaxies in the CMDs redder than the red- 
sequence indicates that the contamination rate is small. 
We find that group galaxies exhibit color-magnitude dis- 
tributions similar to clusters: that is, group galaxies can 
be separated into two distinct populations of red and 
blue galaxies, and the red members, including those in 
poor groups, form a clear sequence. 

We fit the red sequence using red galaxies with 
log(M*/M0)>lO.6, and present the resulting fit in Table 
2. We use a relatively high stellar-mass limit to perform 
the fit to minimize possible confusion arising from the 
blue cloud galaxies due to larger photometric uncertain- 
ties at the faint end. These red group galaxies are defined 
in the same manner as described in §3.7. We further 
impose a red-cut of 0.3 mag redder than the model red- 
sequence color, in order to minimize the effect on fitting 
from the small number of galaxies in the group galaxy 
sample which are are likely background galaxies with col- 
ors significantly redder than the red sequence model. Be- 
cause of the more sparse data for the low stellar-mass 
groups, we fix the fitting slope for these groups to be 
that of the average of the two bins of the more massive 
groups at the same redshift to derive the red-sequence 
colors and dispersion. We find that the slopes and zero 
points of the fitted red sequences are similar among dif- 
ferent redshift and M^ ^^p bins - an indication that the 
formation epochs for M* galaxies in these groups are 
typically at z ?:2 (e.g., Stanford et al. 1998; Gladders et al. 
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TABLE 2 
The Fitted CMR for Fig. 9 
IN THE FORM AS {/ = A + Bx OR {/ = C + B(a; + 21.41) 



z 




^'i* , Q rp 


A 


B 


C 


0.00- 


-0.35 


12.2- 


-14.0 


0.963±0.049 


-0.029±0.002 


1.584 


0.00- 


-0.35 


11.8- 


-12.2 


0.723±0.068 


-0.039±0.003 


1.558 






11.2- 


-11.8 


0.796±0.001 


-0.034* 


1.524 


0.35- 


-0.45 


12.2- 


-14.0 


0.834±0.039 


-0.035±0.002 


1.583 


0.35- 


-0.45 


11.8- 


-12.2 


0.520±0.098 


-0.049±0.004 


1.569 






11.2- 


-11.8 


0.605±0.002 


-0.042* 


1.504 


0.45- 


-0.55 


12.2- 


-14.0 


0.274±0.075 


-0.061±0.003 


1.580 


0.45- 


-0.55 


11.8- 


-12.2 


0.071±0.162 


-0.069±0.007 


1.548 






11.2- 


-11.8 


0.098±0.005 


-0.065* 


1.490 



using the mean slope of the other two Mt.grp bins at the same 
redshift. 



1998). The average slope of ~ —0.046 for the 6 more nias- 
sive group samples is consistent with the value of -0.047 
in the model of Kodama & Arimoto (1997). We see a 
consistent trend of increasing slope, by about -0.03 from 
the z = 0.25 to z = 0.5 bins, which is entirely in agree- 
ment with that observed for clusters over this redshift 
range (but for slightly different filters) in Gladders et al. 
(1998, see their Figure 4) and their models. 

The average red-sequence rest-frame zero point {B — 
Rc)q at M'fj^ = —21.41 is 1.570 mag, and has a small 
root-mean-squared (rms) dispersion of 0.016 mag over 
all redshifts and M+^g^p bins. However, it appears that 
there is a trend of a slightly bluer {B — Rc)o zero point 
for less massive groups: The zero points averaged over 
redshift are 1.571, 1.576, and 1.564 for the massive-, 
intermediate-, and low-mass group samples. These dif- 
ferences account for almost all the dispersion in the red- 
sequence color zero points, as the variations over redshift 
within each group mass grouping are only ~ ±0.01 mag. 
It is not clear whether this small observed difference is a 
reflection of a trend in the average age since formation 
or the average metallicity in the luminous red-sequence 
galaxies in groups of different masses. 

We also compute the dispersion of the rest-frame 
B—Rc colors for the bright (< M*) red-sequence galaxies 
in the 9 subsamples. They range from 0.08 to 0.10 at the 
low- and high-redshift bins for the more massive groups, 
to about 0.10 for the low- mass groups at all redshift bins. 
However, there are two observational uncertainties which 
contribute to the dispersion of the red-sequence. First, 
the photometric redshifts of the groups have a typical 
uncertainty of 0.03, which translates into an uncertainty 
of about 0.03 mag in the k-corrections for red galaxies. 
Second, the typical mean photometric uncertainties in 
B — Rc range from 0.03 mag for the low-redshift samples 
(at the brighter end of our photometry) to 0.06 mag for 
the high-redshift samples. Rough corrections in quadra- 
ture for the dispersion due to these effects give the in- 
trinsic red-sequence color dispersion to be ^ 0.03 for all 
cases. While the rough corrections and the possible con- 
tamination from background galaxies prevent us from 
discerning any pattern in the red-sequence dispersions 
for the different subsamples, the estimated average dis- 
persion is consistent with observations of much smaller 
samples of rich clusters with more precise redshift and 
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Fig. 10. — Red-galaxy fraction, fred, as a function of redshift in 
three log(M»^jrp) bins using all groups galaxies within 0.5i?200 
without and with stellar mass controlled. Groups with larger 
M»,grp exhibit larger /^^j, indicating an apparent group down- 
sizing effect. 

photometric information (e.g., Sandage & Visvanathan 
1978; Bower et al. 1992; Ellis et al. 1997) 

This high degree of uniformity in the red sequence in 
our group samples suggests that the more massive (~ 
M^^, or log(M*/MQ)> 11.0) red-sequence galaxies are 
very similar by z ~ 0.5 in all groups with a halo mass 
larger than approximately IO-'^^'^Mq; i.e., their formation 
time is likely to be at z > 2, similar to the conclusions 
drawn from more massive clusters (e.g., Gladders et al. 
1998; Quadri et al. 2007). From a visual inspection of 
the CMDs, the number of blue group galaxies appear 
to drop significantly from z ~ 0.50 to z ~ 0.25. In 
the next section, we quantitatively measure the fraction 
of red group member galaxies (/red) and examine their 
dependence on group and galaxy properties. 

4.2. Dependence of fred on Group Total Stellar Mass 

We probe the red galaxy fraction fred in galaxy groups, 
computed within 0.5i?2ooi a-s a function of redshift for the 
three different M^^grp samples and present the results in 
Figure 10. Panel a) of Figure 10 shows the change in 
fred as a function of redshift to z = 0.45 using a galaxy 
stellar-mass limit of log(M*/M0)=lO.2. Note that we 
measure fred only out to z 0.45, as at z ~ 0.5 the 
galaxy sample is not complete at the log (M*/Mq)=10.2 
level. All three M^^g^p subsamples exhibit increasing 
fred toward lower redshift; i.e., the Butcher-Oemler ef- 
fect (Butcher & Oemler 1984) observed on the scale of 
galaxy groups. The more massive groups have a larger 
fred in all redshift bins. In other words, more massive 
groups turn their galaxies red, due to the suppression of 
star formation activity, at an earlier time than those in 
less massive groups. This result can be described as a 
down-sizing evolution effect of sorts operating in galaxy 
groups. 

To ascertain that this effect is not the result of rel- 
atively more massive galaxies inhabiting more massive 
groups, we further control for the stellar mass of group 
member galaxies in the next three panels of Figure 10. 
To aid in the discussion of the evolution of fred, we 
use /rT=0.85 as the threshold fred for which a sam- 
ple of galaxies is considered as having the bulk of its 
galaxies with their star formation quenched. Panel b) 
shows that the most massive group member galaxies 
(log(M*/MQ)>ll) exhibit essentially flat fred-z trends 
for groups in different M^ g^p bins with fred>frT, with 
the exception of the z ~ 0.5 point for the low-mass 
groups. That is, these massive group galaxies have 
mostly finished their star formation by z ~ 0.5, regard- 
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less of the influence of group environment. 

In contrast, the effect of M^^grp is clearly observed for 
the subsamples of galaxies of intermediate and low stel- 
lar masses. First, galaxies of similar stellar masses are 
redder in more massive groups (panels c and d), regard- 
less of redshift, demonstrating the dependence on M^^g^p- 
Second, the fred-z trends are different for intermediate- 
and low-mass galaxies in groups of different M^^g^p- The 
intermediate-stellar-mass galaxies in the massive groups 
reach frT at z between 0.45 to 0.25, depending on the 
mass of the groups to which the galaxies belong. For the 
low-stellar-mass galaxies, only those in the most mas- 
sive groups manage to barely reach frT at z ~ 0.25. In 
general, the low-mass galaxies also show larger changes 
in fred than the intermediate-stellar-mass galaxies over 
z ^ 0.45 to z 0.25 range. The different behaviors 
of the fred-z trcuds for these three subsamples indicate 
that, while the fred correlates strongly with the stellar 
mass of the group galaxies, the halo mass of the group in 
which the galaxies are situated also has a significant in- 
fluence on the evolution of the lower-mass group galaxies. 
Galaxies in higher mass groups reach frx at an earlier 
epoch and have a less steep fred-z trend. 

The large fred and flat trends with redshift for the mas- 
sive group galaxies further provides evidence of galaxy 
'down-sizing' studied in the literature; i.e., star forma- 
tion takes place and also stops earlier in massive galaxies, 
then shifts to less massive systems at later times (e.g., 
Cowie et al. 1996; Cattaneo et al. 2008; Seymour et al. 
2008). Massive galaxies, by and large, show little effect 
from being in groups of different masses (at these mod- 
erate redshifts) ; their evolutionary stage is likely a reflec- 
tion primarily of their own mass. However, lower-mass 
galaxies do show evolutionary histories significantly cor- 
related with the halo mass of their parent groups. Hence, 
the 'group down-sizing' observed in panel a) of Figure 
10 is not caused by differences in the stellar mass dis- 
tributions of member galaxies in groups with different 
richness, but rather is the result of lower-mass member 
galaxies in higher-mass groups being in a more advanced 
evolutionary state than those in lower-mass groups. 

4.3. Dependence of fred on Group-Centric Radius 

In the literature, many studies have suggested that 
group environment can truncate star formation in their 
member galaxies (e.g.. Font et al. 2008; Verdugo et al. 
2008; van den Bosch et al. 2008; Wilman et al. 2008, L09). 
For example, Domfnguez et al. (2002) have shown, using 
groups from the 2dF, that the fraction of galaxy popu- 
lations in groups, especially the rich and massive ones, 
exhibit similar radial trends as those observed in galaxy 
clusters with an offset in the values of fred- Using the 
CN0C2 group catalog, Wilman et al. (2005) found that 
the fraction of passive spiral galaxies, which are galaxies 
of spiral morphology but with no/little [OH] emission 
lines, increases continuously toward the group center, 
and demonstrate that star formation rate is truncated 
gradually as galaxies move toward the group center. 

To quantify radial color trends within galaxy groups, 
we compute fred as a function of Vgrp. We separate the 
group sample by group total stellar mass M^^grp, since we 
have found in §4.2 that the fred within groups at a fixed 
redshift depends on M,^grp- The results, using group 
galaxy samples with a stellar mass limit of log(M*/M0)> 
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Fig. 12. — '/red-f'grp 'trends for group galaxies with log(M»/Af0) 
> 11.0 (solid pink), log{M,/M0)=lO.6-ll.O (dotted purple), 
and log(Mt/MQ)=10.2-10.6 (dashed blue) in each redshift and 
log(M*,9rp) bin as indicated in the panels. A simple linear fit 
is applied to each trend with the slope indicated. In general more 
massive group member galaxies exhibit gentler slopes toward to 
the center, and more massive groups produce steeper radial depen- 
dence. 

10.2 are presented in Figure 11 for the two lower redshift 
bins. 

Generally speaking, group galaxies in all M^^grp and 
redshift bins exhibit declining fred with increasing rgrp. 
The central regions of the groups are dominated by red 
members with fred^ 0.6 to 0.9, dropping to '--^0.3-0.4 at 
the outskirts. For both redshift bins, the effect of M, 
is clearly observed within rgrp sl, where groups of higher 
^*,grp have larger fred- Beyond i?200; the differences 
in fred among different M»^grp groups become smaller, 
and approaching similar values of fred'^ 0.3 to 0.4 for all 
groups. However, the observed trends in Figure 11 are 
likely driven differently by group galaxies of different stel- 
lar masses. We accordingly re-plot the fred-Tgrp trends 
by dividing our group galaxies into three subsamples, 
with log(M,/Mo) > 11.0, log(M,/M0) = 10.6 - 11.0, 
and log(M,/M0) — 10.2 — 10.6. For the two subsamples 
of more massive galaxies, we also include the z ~ 0.5 
redshift bins. The results are presented in Figure 12. As 
an aid to interpretation, we also perform linear fits of 
the fred-fgrp trcuds for the 24 subsamples, with the best 
fitting slopes shown in the Figure. 

While the error bars are large for these plots, sev- 
eral trends are apparent. As expected, the most massive 
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Fig. 13. — Local galaxy density S5 as a function of rgrp in groups 
with different M*^g,.p. Larger circle sizes indicate galaxies at lower 
redshift. The solid curves are the fitted 2D NFW (Navarro et al. 
1996) relations using the analytic form in Bartelmann et al. (1996). 

galaxies (log(M»/MQ)>ll) show little effect from their 
location within their parent group, as almost all have es- 
sentially fred^frT- The changcs in the fred-Tgrp trends 
with z or M*_grp are primarily driven by group galax- 
ies with intermediate (10. 6< log(M*/MQ)<11.0), and 
low masses (10.2<log(M*/MQ) <10.6). These fred-fgrp 
trends follow a couple of general patterns. First, the low- 
mass galaxy samples have the steepest radial changes of 
fred for all group mass and redshift bins, indicating that 
they are affected by the group environment the most. 
Second, the radial slope of fred appears somewhat shal- 
lower for these galaxies in low-massive groups, which is 
suggestive that more massive halos may have a relatively 
larger effect on the galaxies that they are accreting, al- 
though larger samples are needed to verify this effect. 

4.4. Local Galaxy Density for Group Galaxies 

It has generally been accepted that local galaxy den- 
sity can affect galaxy properties such as star forma- 
tion rates and colors (see §1). Using a small sample of 
rich clusters, L09 probed the local galaxy density effect 
on cluster galaxies at different cluster-centric radii over 
0.15 < z < 0.55. They showed that local galaxy density 
has an important effect at cluster outskirts but less so 
in the cluster central regions. Here, with a much larger 
sample covering a larger range of halo masses, we exam- 
ine the interplay between the local galaxy density effect 
and group environmental influence, as parameterized by 
the group-centric radius, on group galaxies. These two 
measures can be considered different in that the former 
tells us how close neighbors affect the galaxy population, 
while the latter delineates the effect due to the position 
of the galaxy in the larger parent halo. 

Figure 13 plots the local galaxy density, Ss, for each 
group galaxy, which is measured based on the nearest 
fifth neighbor distance using the method in L09, as a 
function of the group-centric radius for groups with dif- 
ferent M^^grp- The distributions cover over three orders 
of magnitude in local density and can be described in 
general as having a decreasing average with increasing 
radius. In each M^^g^p bin, we also show that the lo- 
cal density distribution roughly follows a projected NFW 
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Fig. 14. — Red-galaxy fraction fred within 1.5i?200 a-s a function 
of redshift for galaxies in different local galaxy density S5 subsam- 
ples. The 9 panels separate subsamples of different group stellar 
masses yit,grp and galaxy stellar masses M,, as indicated on each 
panel. The red solid, gray dotted, and blue dashed curves represent 
high-, intermediate-, and I0W-S5 divisions, respectively. 

profile (Navarro et al. 1996), fitted following the prescrip- 
tion of Bartelmann et al. (1996). These plots show that 
the E5 measurement produces meaningful relatively local 
density information when the density range is sufficiently 
large, e.g., over a factor of several from low to high den- 
sity. Typically, the group core has local galaxy densities 
that are close to 100 to several hundred times higher 
than those at 2i?2oo, for low- and high- mass groups, re- 
spectively. Within 0.5i?2oo, S5 decreases smoothly with 
Tgrp. Beyond that, some regions with higher E5 are 
super-imposed on the envelope, indicating possible in- 
falling sub-groups at large Vgrp^ or the ambiguity in group 
membership among multiple groups. 

4.4.1. The Effect of Local Density 

To examine the effect of local galaxy density environ- 
ment on the galaxy population, in Figure 14 we plot 
fred as a function of redshift for group galaxies within 
1.5i?2oo into three bins of local galaxy density. We fur- 
ther separate the galaxy samples into bins of galaxy stel- 
lar mass and total group stellar mass as in §4.3, since 
we have found that fred has a dependence on both of 
these parameters. Figure 14 provides a look at the ef- 
fect of E5 integrated over the whole group. Our highest 
and lowest S5 bins contain the highest 20-percentile and 
the lowest 40 percentile of group galaxies in our sample. 
The median local galaxy density values of the high- and 
intermediate-density bins are approximately 12.8 and 3.7 
times that of the low-density bin, respectively. 

The main feature of Figure 14 is that the effect of 
S5 depends on galaxy stellar mass. Within the factor 
of ~13 between the low- and high-density bins, fred in- 
creases by ~ 0.2 for the intermediate-mass galaxies, and 
~ 0.4 for the low-mass galaxies; whereas for the high- 
mass galaxies, there is no statistically significant differ- 
ence in fred- The trends for the high- and intermediate- 
mass galaxies appear to be very similar for galaxies in 
groups of different M*^grp, and also over the redshift 
range of 0.5 to 0.25. For the low-mass galaxies, the dif- 
ferentials in fred due to different E5 may be somewhat 
larger for galaxies in massive groups. Thus, the mag- 
nitude of the local galaxy density effect appears to be 
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Fig. 15. — Red-galaxy fraction /j-e^ a-s a function of group-centric 
radius Vgrp for group galaxy subsamples of different local galaxy 
density E5. The subsamples are further divided by redshift and 
group member galaxy stellar mass M* bins, as indicated on each 
of the 8 panels. The different S5 subsamples from high to low 
are plotted in different linestyles as (red) solid, (gray) dotted, and 
(blue) long-dashed, respectively. Groups of all total group mass 
M*,grp are combined for better signals. 

fairly similar for galaxies in groups of different richness 
across the redshift range of 0.5 to 0.2, but has a strong 
dependence on the stellar mass of the galaxies. 

4.4.2. Local Galaxy Density and Group Environment 

The parameters E5 and rgrp are correlated. To exam- 
ine how these two effects interplay with each other, in 
Figure 15 we present fred of samples of group galaxies of 
different stellar masses as a function of rgrp for galaxies 
in three E5 environments at different redshifts. Since it 
appears that the group halo mass has only a weak influ- 
ence on the effect of S5 on fred, we combine groups of 
different masses in order cut down the uncertainties in 
the fred measurements, which are still fairly large when 
the galaxy samples are further divided up into rgrp bins. 

Inspection of Figure 15 shows that at all group-centric 
radii, the local galaxy density has the strongest influence 
on the low-mass group galaxies, as is the case in Figure 
13. In general, the effect of E5 on fred is similar for 
galaxies at different group-centric radii. The most mas- 
sive galaxies subsamples, again, all have high fred (larger 
than 0.8) for all local densities and group-centric radii. 
The fred of intermediate-mass galaxies shows a possible 
weak change with local galaxy density with low signifi- 
cance, but consistent with that measured in Figure 13. 
For the low-stellar mass group galaxies there is a clear 
correlation between fred and local galaxy density. This 
is observed in all group-centric radii, including the core 
region. 

Importantly, note that for the low-mass galaxies, even 
with local galaxy density controlled, galaxies are on av- 
erage redder in the center of the group than those in the 
outskirts. A similar, but weaker, trend is also seen for the 
intermediate-mass galaxies. This strongly suggests that 
the location of the (lower-mass) galaxies in their parent 
group halo has an effect on their evolutionary history, 
after the effects of local density and stellar mass are con- 
trolled. That is, the redder population seen in the center 
of galaxy groups and clusters is not due entirely to these 



galaxies being in a denser region; being in the center of 
a massive halo also contributes to them turning red. 

5. DISCUSSION 

5.1. The Uncertainties 

Our work is exploited by photo-z, which is known to 
have less accurate redshift information than spectroscop- 
ical measurement. With intricate corrections in weight 
and background subtraction, we discuss uncertainties in 
our procedure and possible effects on our results with 
the aid of mock catalog whose properties are detailed 
in Croton et al. (2006). We assign photo-z to each mock 
galaxy such that the overall photo-z dispersion compared 
to mock redshift is comparable to that in our training set. 

5.1.1. Wi 

For each galaxy whose photo-z uncertainty is less than 
0.6(1-1-2;), we assign a weight Wi based on the ratio of 
the number of all galaxies in the catalog to that in the 
sample in a given magnitude bin (§2.2). Such Wi is only 
a function of magnitude, because we find that both blue 
and red galaxies, roughly separated hy B — Rc=1.8, have 
comparable Wi . The role of Wi is to account for the miss- 
ing galaxies in a given magnitude bin due to large photo-z 
errors. In practice, for n galaxies drawn from the sample, 
the expected total galaxy count is X]i=o n '^i ■ 

To estimate the effect of Wi , we apply the same Wi com- 
putation to the mock catalog with the simulated photo-z. 
The Wi reaches 2 at i?c=22.45, typical to our observed 
data. We then compute fred for the mock galaxies at 
each redshift bin. We find that such fred is smaller by 
'--^10% than that using the known mock redshifts (i.e., no 
galaxy is removed hence no Wi correction) at each red- 
shift bin. In other words, the correction results in a 
smaller fred but such effect is uniform over different red- 
shift bins. The zeropoint of the fred — z trend is smaller 
with Wi applied, but the slope of the trend remains the 
same. 

5.1.2. Missing Groups and Contaminated pFoF Groups 

Our group sample is constructed using the pFoF al- 
gorithm on the photo-z catalog. No group finding al- 
gorithms are known to be perfect, including the pFoF 
algorithm. We apply the pFoF algorithm to the mock 
catalog with the simulated photo-z and computed Wi. 
We test the performance of the algorithm and investigate 
the roles of missing mock groups and falsely detected and 
contaminated FoF groups. 

First, we find that the pFoF algorithm is able to re- 
cover >80% of the mock groups whose halo mass are 
> 2 X IO^^Mq within our redshift range, and the recov- 
ery rate is not a strong function of redshift. We compute 
fred using mock galaxies brighter than A/^ + 1.5 in each 
mock group, and investigate the average fred with and 
without any missing mock groups. Since the recovery 
rate depends on halo mass, we carry out the compari- 
son in three halo mass bins. We find that for groups 
with halo mass greater than 10^^ Mq, approximately the 
lower limit of our sample, the missing groups have essen- 
tially no effect on the average fred, with the value being 1 
to 2 % redder for the pFoF recovered groups when com- 
pared to the complete mock group catalog. It is only 
by going down to the halo mass range in halo mass to 1 
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to TxlO^^M©, well below our nominal group mass limit, 
where we find an effect of about 10%. 

Second, we find the fraction of false and contaminated 
mock pFoF groups is a strong function of redshift and 
M^^grp. We define a contaminated group as a pFoF group 
with a match from the mock group catalog, but with a 
significant number of linked galaxies not belonging to the 
true mock group (see Li & Yee for a precise definition). 
In principle, in a large sample the effect of groups deemed 
contaminated on the measurement of fred is expected to 
be statistically corrected by the background corrections. 
However, groups contaminated with larger fore- or back- 
ground galaxies would be more likely to be identified by 
any FoF algorithm. Thus, we inchide both false detected 
groups and contaminated groups in our tests on their 
effects on the measurements of fred- Hereafter, we will 
simply refer to these groups as 'contaminated groups.' 
The contaminated groups tend to occur at higher redshift 
or small Af^.grp regime. Since galaxies in contaminated 
groups are expected to be bluer on average, contaminated 
pFoF groups makes the fred-z trend slightly steeper in 
the low M^^grp bin. For the low-mass pFoF mock groups 
with M^_grp between 10^^-^ and lO^^-^M©, the contam- 
inated group fraction is ^-^20% at zsO.45, but increases 
to ^ 67% at z ^0.5. Using the mock groups we find fred 
for the mock pFoF sample at z ~0.25 to be a factor of 
^ 1.07 smaller than the true value, with the discrepancy 
increasing to ~ 1.23 smaller at z ~0.5. For the more 
massive pFoF mock groups, the contaminated fraction is 
only ~15% at z ^0.5, resulting in a smaller average fred 
by ^-^3%. We note that this small effect comes primarilj^ 
from contaminated groups, rather than false detection 
groups, as there are negligible number of false massive 
groups. The contamination fractions at lower redshift 
are smaller, resulting in negligible effects on fred for the 
more massive groups. 

We note that these pFoF mock groups are selected by 
Ngai > 5, Ngz > 5 and M^^grp, but not by Bgc. Based 
on Fig. 8, the Bge > 125 cutoff removes 10% of the pFoF 
groups from the sample. Therefore, if they are further 
selected by Bgc, the above false detection fraction and 
the difference in the average fred using the mock pFoF 
groups are expected to be even slightly lower. 

As an example, In Figure 10 panel a), the contam- 
inated pFoF groups in our sample make the fred — z 
trends steeper than the actual slopes. However, based 
on the above mock tests, we argue that the general de- 
clining fred — z trend cannot be simply caused by the 
increase in contaminated pFoF groups, since the differ- 
ence in fred between z ^0.25 and z ~0.5 is considerably 
larger than what is expected due to the effect of contam- 
inated groups. To account for the possible systematic 
effect of false group detection and contaminated groups, 
we increase the upper uncertainty error bars on fred for 
the low-mass group bins based on our mock group tests 
in all Figures involving fred- 

5.1.3. Group Centers and Group-centric Radius 

In §3.3 we assign a pFoF linked galaxy to be the cen- 
ter of that galaxy group. However, it is not necessary to 
have a galaxy located at the center. We do so because 
this allows us to study properties of these 'central' galax- 
ies in the future. Another way to define the group center 
is to take the median position of all the linked galax- 



ies weighted by their luminosity and local galaxy surface 
density. L09 find that such centers for the CNOCl clus- 
ters are about 30" from the cD galaxies. For the RCSl 
groups in this work, the median separation of these two 
types of group centers is ^8", which is equivalent to 0.04 
Mpc at z=0.3 or O.OSMpc at z=0.5. 

The i?200 is another factor which affects the group- 
centric radius. We compute the -R200 by extrapolating 
the Bge - i?200 relation in Yee & EUingson (2003). The 
derived i?2oo for our group sample is ~1.07±0.24Mpc. 
Therefore, the uncertainty in ii200 dominates over the 
role of group center in computing Vgrp. We note that 
the Bgc — -R200 relation, which is derived based on the 
CNOCl clusters, may not be the same for galaxy groups. 
However, the relation indeed offers us an efficient way in 
estimating ii200 to be within an order of magnitude to 
the actual values. 

The uncertainty in the i?2oo introduces 
Alog{M^,^grp) ^ O.lA'/0 to the M^,grp in §3.4, and 
gives a Afred ~ 0.04 in the fred computed using group 
galaxies within 0.5i?2O0 of the groups. We believe that 
our conclusions remain the same even though we have 
large uncertainties in group-centric radii, because our 
binsize for M^^g^p is large enough to accommodate this 
extra uncertainty due to i?200) and any changes in fred 
are bigger than 0.04 in general. 

5.2. Properties of Group Galaxies 

It has been known for decades that galaxy properties 
correlate with their environment: Galaxies in clusters are 
dominated by a red population, and field galaxies are 
characterized by a blue population. The morphology- 
density relation presented in Dressier (1980) has set the 
stage for many studies in the following decades to quan- 
tify this relationship and to investigate the causes. With 
the halo model of galaxies and groups, Weinmann et al. 
(2006) suggest that the morphology-density relation can 
be expressed in terms of halo mass instead of galaxy num- 
ber density, because the projected galaxy number in a 
halo is expected to correlate with the halo mass. A stan- 
dard picture of the halo model is that some galaxies are 
embedded in their own individual halos independently, 
and some reside within a common halo shared with other 
galaxies in the form of galaxy groups. Once a halo enters 
within the virial radius of a more massive one, such as 
when a single-galaxy halo enters a group halo, it is re- 
ferred to as a satellite halo or sub-halo within the larger 
halo (or, the parent halo). As the satellite or sub-halo 
orbits within the more massive one, it may be subjected 
to all kinds of environmental effects; its mass is reduced 
and the diffuse outer part can be stripped off by tidal 
effects and interactions with other substructures. In ad- 
dition to the mass loss, the gas reservoir of the infalling 
galaxy will be deprived and eventually lead to a halt in 
its star formation when the galaxy depletes its cold gas. 
This is usually referred to as 'strangulation' (e.g., Balogh 
et al. 1999, 2000; van den Bosch et al. 2008; Kawata & 
Mulchaey 2008; McCee et al. 2009). In this section we 
discuss the results from our galaxy group sample in the 
context of the halo model and 'nature' versus 'nurture' 
scenarios. 

5.2.1. Four Factors that Affect the Group Galaxy Population 
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We have identified one intrinsic and three global pa- 
rameters that affect the evolution of the population of 
group galaxies as measured by the red galaxy fraction, 
fred- The underlying predominant parameter is the mass 
of a galaxy's own halo (measured via the galaxy's stel- 
lar mass). The other three, (1) the mass of the parent 
group halo into which the galaxy has fallen, (2) the po- 
sition of the galaxy in the parent halo, and (3) the local 
galaxy density, produce further effects which evidently 
accelerate the evolution of the galaxies to their final red 
color. Here we summarize the observed effects of these 
four parameters on fred as shown in Figures 10- 15. 

• The galaxy population as measured by fr^d is strongly 
dependent on the galaxy stellar mass; i.e., its own halo 
mass. The most massive galaxies (log(M*/M0)>ll, or 
Mftc for red galaxies) are already mostly red 
in all group environments by z ^ 0.5. The effects of 
the group environment and the local galaxy density, 
summarized below, are primarily seen in galaxies with 
log(M*/M0)<ll. 

• The mass of the parent group halo affects the galaxy 
population such that the red galaxy fraction is larger for 
galaxies in more massive groups. This effect is observed 
after controlling for the stellar mass of the galaxies and 
it is stronger for lower-mass galaxies. 

• The red galaxy fraction has a dependence on the loca- 
tion of galaxies in their parent halo, as parametrized by 
the group-centric radius, rgrp, normalized by i?200- The 
fred for intermediate- and low-mass galaxies increases 
toward the group center, with the radial gradient being 
steeper for lower-mass galaxies. For low-mass galaxies, 
there are also some indications that the fred gradient is 
steeper for galaxies in more massive groups. 

• There is a degeneracy in the two parameters, group- 
centric radius (rgrp) and local galaxy density (S5), that 
measure the environment of group galaxies. By examin- 
ing galaxy samples of intermediate and low stellar masses 
in bins of local galaxy density, we find that: 1) for galax- 
ies in the same local galaxy density bin, there is still a 
radial gradient in fred', and 2) conversely, for galaxies 
at the same group-centric radius, those in higher local 
galaxy density regions have larger fred- These results in- 
dicate that both the position of the galaxy in the parent 
halo and the local galaxy density around each individual 
galaxy have an effect on its evolutionary history. 

In the following subsections, we discuss in more detail 
the dependence of fred on these various parameters, and 
how it fits into the general context of galaxy evolution. 

5.2.2. The Dependence on Galaxy Stellar Mass 

The apparent strength of the dependence of the fred on 
environmental effects, whether it is the mass of the group 

halo, the position the galaxies in the halo, or the local 
galaxy density, is a strong function of the stellar mass of 
the galaxies. The biggest effects are seen in the lowest 
mass galaxies; while for the most massive galaxies, those 
with log(A'/,/Af0) <;llA'/0, little discernible environmen- 
tal effect is found. Furthermore, for these most massive 
galaxies the fred, which is close to 1, has no apparent de- 
pendence on redshift (within our redshift range of z sO.5) 
regardless of their group environment. The uniformity of 
the slope, color zero-point, and dispersion of the stacked 
red-sequence (Figure 9) for galaxies brighter than 
also points to the conclusion that the most massive galax- 



ies in our sample have completed their evolution to a 
quiescent state by z ^ 0.5, suggesting that they mostly 
formed early and have little or no star formation since 
2; ~ 1. Our conclusion that the stellar mass of a galaxy is 
a more fundamental factor in galaxy evolution echoes the 
recent spectroscopic results (e.g., Griitzbauch et al. 2010; 
Cucciati et al. 2010; Li et al. 2010; Sobral et al. 2010). 
This illustrates that our photometric-redshift samples are 
able to give reliable conclusions while it provides a group 
sample size much larger than any from spectroscopic sur- 
veys at 2; > 0.15. 

However, that the lowest-mass galaxies appear to be 
most affected by the group and local density environ- 
ments docs not necessarily mean that these; environmcm- 
tal factors produce disproportionately larger effects on 
low-mass galaxies, or that massive galaxies are immune 
to them. Rather, it is likely that massive galaxies are 
equally affected by their environment, but the (secular) 
galaxy down-sizing evolution effect has already turned 
most massive galaxies red well before z ~ 0.5. And once 
they are red, the 'nurturing' influence can no longer be 
discerned-whether a galaxy is far away from a parent 
halo, or in a low galaxy density region, it is already red. 
The same environmental effects likely have operated on 
these massive galaxies earlier on, accelerating the process 
of them turning red, to different degrees for galaxies in 
different environments. In this interpretation, we should 
be able to observe similar kinds of dependence of fred on 
group-centric radius, group halo mass, and local galaxy 
density for galaxies of larger masses at higher redshifts. 
By the same token, for intermediate-mass galaxies at 
lower redshifts, we should expect them to have a smaller 
apparent dependence of their galaxy colors on group or 
local galaxy density environments, because many of them 
would have already turned red by the present epoch. 

5.2.3. The Dependence on Group Halo Mass 

In the halo model, the infall of a galaxy in a subhalo 
into the parent group halo is a significant event in the 
evolution of the stellar component of the galaxy. The 
larger halo can affect the star formation history of the 
infalling galaxy via a number of mechanisms that act on 
a global scale within the parent group, such as strangu- 
lation, ram-pressure stripping, and tidal disruption. The 
dependence of the galaxy population on the mass of their 
parent group or cluster can provide insights into the ef- 
ficiency of the various environmental mechanisms that 
affect star formation history, and their relative impor- 
tance. There is some evidence from redshift surveys in 
the local universe that galaxy population properties have 
a dependence on the mass of the parent group halo, in 
that galaxies in more massive halos tend to be in a more 
evolved, or quiescent, state (e.g., Weinmann et al. 2006; 
Kimni et al. 2009; van den Bosch et al. 2008, L09). At 
higher redshifts, lovino et al. (2010), using groups with 
z up to 1 cataloged from the zCOSMOS survey, found 
that the galaxy blue fraction is larger for groups of lower 
richness. We find in our sample a similar dependence in 
general, and we examine this dependence in more detail 
as functions of both galaxy stellar mass and group-centric 
radius. 

The dependence of the galaxy population on the parent 
group halo mass is primarily seen in lower-mass galaxies. 
Figure 10 also shows that the effect has a redshift de- 
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pendence. By z ^ 0.25, the difference in /^ed for galax- 
ies in groups of different masses is seen mostly in low- 
mass galaxies; whereas at z > 0.35, both intermediate- 
and low-mass galaxies show a significant dependence of 
fred on M*_grp- Further subdividing the samples by the 
group-centric radius, we find that the increase in j^ed for 
the intermediate- and low-mass galaxies occurs over the 
range of group-centric radius from the core to outside the 
virial radius. However, in the core of rich groups galaxies 
of different masses have a narrower range of fred values, 
compared to those in poor groups; i.e., the dependence 
on group-centric radius is stronger for galaxies in rich 
groups than poor groups. We note that some caution 
should be \ised when considering the strength of this ef- 
fect, in that it is more difficult to define the center of a 
poor group, which may cause a smearing of the radial 
dependence. Weinmann et al. (2006) observed a simi- 
lar dependence of the fred-fgrp trend on group mass at 
z ~ using their SDSS group sample, with their lowest 
mass groups (which are less massive than our low group 
mass bin) showing almost no dependence of fred on rgrp. 
Over all, this dependence on the group halo mass can 
be described, to borrow a popular term from galaxy evo- 
lution, as some sort of 'down-sizing' effect, in that the 
galaxy population in richer groups appears to be on av- 
erage in a more advanced stage of their evolution. 

Our data demonstrate the additional influence that 
group environment has on their constituent galaxies, but 
does the effectiveness of this influence depends on the 
mass of the parent halo? That is to say, is this apparent 
'group down-sizing' effect due to the different effective- 
ness (or prevalence) of the mechanisms that quench the 
star formation of the accreted member galaxies? While 
this may be the case, the same effect can also be pro- 
duced by a simpler explanation: The larger fred in more 
massive groups can simply be a result of a hierarchical 
build-up of structures, such that at any given epoch the 
more massive groups/clusters have been in existence for a 
longer time than less massive groups. In a scenario where 
galaxy stellar mass is the primary determinant of the 
time scale of galaxy evolution, we can consider that group 
environment is acting as an accelerator for the galaxy to 
reach its final red quiescent state (e.g., Bolzonella et al. 
2010; Balogh et al. 2011). At a fixed redshift, ff the time 
scale for the growth of the group from the accretion of 
'field' galaxies is longer than that for the truncation of 
star formation in the galaxies, then the infalling galaxies 
would turn red faster than the rate of accretion of more 
field galaxies, creating a situation where groups that have 
been in existence for a longer time would have a larger 
fraction of red galaxies. Under this scenario, clusters 
or groups of a fixed halo mass would expect to have an 
overall bluer population at higher redshift than those at 
lower redshift, since, even if their accretion histories are 
similar, their accretion time scale is compressed (McGee 
et al. 2009) . This situation would naturally produce the 
Butcher-Oemler effect for a sample of fixed group/cluster 
mass over a range of redshift. 

Observationally it is impossible to trace precisely the 
evolution of the progenitors of a given group along 
the redshift axis. Hierarchically, structures grow from 
smaller systems into larger ones; a rich group at lower 
redshift is likely to have a smaller mass in the past. In 
Figure 10, we use fixed cutoffs at different redshift 



bins without considering such mass evolution. If we could 
trace the progenitors, i.e., distinguish the ancestors, the 
apparent Butcher-Oemler effect and the 'group down- 
sizing' effect is expected to be even stronger. For exam- 
ple, a group at 2; ~ 0.25 is expected to have a smaller 
mass at 2; ~ 0.5, and hence a lower M^^grp. Since groups 
of smaller M*,gj.p at a fixed redshift have smaller fred, the 
'fred-z' gradient for this group in our example will be- 
come even steeper, making the 'group down-sizing' effect 
more significant. 

5.2.4. Environmental Influences: The Effect of Local Galaxy 
Density and Group-Centric Radius 

The local galaxy density (S5) and group-centric (or 
cluster-centric) radius (vgrp) are two parameters that are 
often used to measure the environment of galaxies. In 
fact there has long been a running discussion in the lit- 
erature as to which is more fundamental in determining 
the galaxy population properties of cluster galaxies (e.g.. 
Dressier 1980; Whitmorc et al. 1993). While they are 
correlated with each other, since local galaxy density is 
higher in the center relative to the outskirts, these two 
parameters can provide indications of different physical 
mechanisms. The group-centric radius can be consid- 
ered as a parameter sensitive to the global environment 
of the gravitational potential and gas content of the par- 
ent group halo. Thus, Vgrp can be used to gauge the 
importance of mechanisms such as ram-pressure strip- 
ing, strangulation, and global tidal effects from the group 
dark matter halo mass. On the other hand, local galaxy 
density measures the more immediate environment, and 
is more pertinent to mechanisms that are dependent on 
the existence of neighboring galaxies, such as harassment 
and galaxy-galaxy interactions and mergers. 

In most studies of group or cluster galaxies, where the 
galaxy samples arc usually too small to separate into 
subsamples of both radial and local-density bins, corre- 
lations of galaxy population properties with one or the 
other parameter become ambiguous. Is the change in the 
galaxy population as a function of group-centric radius a 
reflection of the different local galaxy densities at differ- 
ent radii? Or conversely, is the correlation between local 
galaxy density and galaxy population primarily a result 
of the position of the galaxies within the larger parent 
halo in which they reside? These questions can be an- 
swered by examining the properties of galaxies at a fixed 
group-centric radius as a function of local galaxy density, 
or for galaxies in regions of a specific local galaxy den- 
sity in different group-centric radii. A similar but more 
bimodal view is taken by some investigators in which 
galaxies are simply divided into "central" galaxies and 
"satellite" galaxies (e.g., van den Bosch et al. 2008). 

Figure 15, which plots the fred of group galaxies in 
different local galaxy density bins as a function of group- 
centric radius, demonstrates the residual effect of local 
galaxy density after accounting for the stellar mass of the 
group galaxies and their positions within their parent 
groups. Again, the low-mass galaxies show the largest 
effect (showing up as the separations between points of 
the same symbols in each panel of Figure 14) due to 
local galaxy density, where over the factor of about 13 
in S5, fred is changed by roughly 0.2 to 0.5, with the 
difference being more significant at lower redshift. This 
change is comparable to the changes seen in fred over the 
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gc from outside i?2oo to tire core. Similarly, for 
the intermediate- and high-mass galaxies, the difference 
in fred induced by the local galaxy density is also the 
same order as that seen for these galaxies between ~i?200 
and the group core: A fred < 0.1 for the massive galaxies, 
and about 0.1 to 0.2 for the intermediate galaxies. 

In most stiidies of the effect of local galaxy density on 
the make-up of galaxy populations, the location of the 
galaxy within a larger halo is rarely considered (since of- 
ten, these are "field" samples). Furthermore, the stellar 
mass effect is clearly dominant over that of local galaxy 
density. Thus, if neither of these are controlled, the lo- 
cal galaxy density correlation could become misleadingly 
strong. There are several investigations of the effect of lo- 
cal galaxy densities for field galaxies in which the mass or 
luminosity of the galaxies is taken into consideration. Us- 
ing the same RCSl photometric redshift catalog as a field 
sample between z of 0.2 and 0.6, Yee et al. (2005) found 
almost no dependence; of /, (;,; on local galaxy density at 
z > 0.5, and a weak dependcince at z ~ 0.3, whereas 
Balogh et al. (2004) found using the SDSS at z ~ a 
strong local galaxy density dependence of fred for galax- 
ies of all luminosities. More recently, Baldry et al. (2006) 
and Peng (2010), using SDSS data, show that most of the 
dependence of fred on environment primarily occurs for 
galaxies with log(M,/MQ)< 10.7. 

Our results, after controlling for group environment, 
show a similar trend to that found by Baldry et al. (2006) 
for SDSS galaxies. Because none of the three parame- 
ters, M*, S5, and fred, are computed identically in dif- 
ferent studies, we can only compare the dependence of 
fred from the two studies approximately. Our three den- 
sity bins center cover a range of approximately 1 . 1 dex in 
S5, each with a half bin width of approximately 0.27 dex. 
It is reasonable to match the local galaxy density param- 
eters, (S5 in this paper, and S in Baldry et al. (2006)), 
by assuming that the high ends of the density measure- 
ments are similar. Hence, the high- and low-density bins 
in our work correspond to approximately logS ~ 0.8, and 
-0.3 in Baldry et al. ^Prom their Figure 11, and aver- 
aged over similar bins of galaxy stellar mass, the change 
in fred, A fred, between local density bins equivalent to 
our low and high bins are ^ 0.1, 0.2, and 0.3 for stellar 
masses of log(M,/M0)>ll.O, 10.6<log (M,/M0)<11.O, 
and 10.2<log (M*/M0)<1O.6, respectively. These rela- 
tive changes in fred for different local galaxy densities 
from the SDSS field galaxy sample are consistent with 
what we find for the group galaxy sample at z 0.25, 
with Afred ^0.1, 0.25, and 0.4 for the three stellar mass 
bins, respectively. 

We also find that the relative effect of local density on 
fred is roughly constant as a function of group-centric 
radius and the mass of the parent group halo. This, 
along with the general agreement with the effect for field 
galaxies, suggests that we can separate the influence be- 
tween local galaxy density and the global effects of the 
group environment. In the context of the debate of 
whether galaxy population is correlated primarily with 
local galaxy density or cluster/group-centric radius, our 
results indicate that both parameters are important. 

5.2.5. Galaxy Groups as the 'Nurturing' Influence 

In early discussions of Dressler's work on the 
'morphology-density' relation (Dressier 1980), along with 



the discovery of the Butcher-Oemler effect (Butcher & 
Oelmer, 1978), attempts to interpret the results along 
the lines of the idea of 'nature' versus 'nurture' have of- 
ten been used, though not much has been resolved (e.g., 
Poggianti et al. 1999; Cooper et al. 2007; Mateus et al. 
2007; Tasca et al. 2009). ^From the perspective of the 
halo model, we can invoke this classic framework in the 
following way. The strong dependence of a galaxy's star 
formation history on its own halo mass can be attributed 
to 'nature;' while how it interacts with its immediate 
neighboring subhalos and how it is being influenced by 
the larger group halo into which it has been accreted can 
be thought of as the 'nurturing' processes. (There can be 
some debate as to whether the local galaxy density can 
be considered as part of the 'nature' process; it could 
be that a galaxy is more likely to be massive because it 
was formed in a high-density region.) The dependence 
of the red galaxy fraction on the stellar mass of galaxies 
can be seen as a manifestation of the currently acccipted 
down-sizing star formation paradigm, and considered as 
the underlying secular evolution of a galaxy. 

After controlling for both galaxy stellar mass and local 
galaxy density, our data clearly show a residual depen- 
dence of the red galaxy fraction on the group environ- 
ment, especially for lower-mass galaxies. Group galaxies 
show a fred gradient with the group-centric radius, in 
that galaxies in the center are redder on average than 
those further out. Furthermore, this effect is larger in 
more massive groups. This gradient is not entirely due 
to the fact that local galaxy density is higher in the core 
of a group, since it is still seen using galaxy samples in 
fixed bins of S5. The change in fred going from out- 
side R200 to the core is comparable to or larger than the 
changes seen in fred due to a change of local density of a 
factor of ~ 13 times. This suggests that the efficiency in 
truncating star formation from physical processes arising 
from close neighbors is of the same order as that due to 
the processes operating in the larger halo into which the 
galaxy has been subsumed. 

Physical processes that are related to neighboring 
galaxies include galaxy mergers and harassment, produc- 
ing strong tidal effects on the galaxy. We can consider 
the local galaxy density effect to be part of the secular 
evolution, on top, or part, of the down-sizing effect; i.e., 
these effects occur independently of whether a galaxy 
with certain local galaxy density is sitting in the field, or 
in the outskirts, or within the virialized region of a larger 
groiip. However, the much more massive group halo may 
produce ram-pressure stripping due to intra-group gas, 
or a strangulation effect from the removal of the outer 
gaseous halo of the infalling galaxy. Such processes can 
serve as additional mechanisms in driving the galaxy evo- 
lution towards the red-and-dead state. The net effect is 
an acceleration of the truncation of star formation as the 
galaxy settles into the parent halo, turning it red earlier 
than field galaxies of similar stellar masses which are not 
affected by a group halo. 

In the current literature, a number of investigators 
have separated the effect of stellar mass and environ- 
ment. For example, Peng (2010) modeled the evolu- 
tion of galaxy population using the concept of mutu- 
ally independent "mass-quenching" and "environment- 
quenching" efficiencies. McGee et al. (2009) used the 
idea of "fraction of environmentally affected galaxies", 
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based on accretion liistory of galaxy groups and clusters 
from n-body simulations, to model the evolution of clus- 
ter galaxy populations. We note that these two studies 
effectively focused on two different environmental param- 
eters - Peng et al. on local galaxy density, and McGcc 
et al. on group environment. In the analysis of our large 
sample of groiip galaxies, we find that we can in fact sep- 
arate the environmental effect into a local galaxy density 
effect and a more global large scale structure effect due 
to the presence of a larger group mass halo into which 
the subhalo has been accreted. These two types of ef- 
fects likely invoke different mechanisms that lead to en- 
vironmental quenching of star formation in the accreted 
galaxies, and both nc^c;d to be accounted for explicitly in 
modeling galaxy evolution. 

In a picture where the the dependence of star forma- 
tion rate on stellar mass is considered as the long term 
secular trend in the evolution of a galaxy, the local galaxy 
density environment can accelerate the quenching rate of 
star formation, but over different time scales, depending 
on the process. For example, galaxy mergers can act as 
a one-time-only process that happens stochastically and 
quickly change the nature of the galaxy; whereas galaxy 
harassment can be a more gentle process acting over a 
longer time scale. The other major environmental event 
is the accretion of a galaxy into a larger halo, which 
will further act as an accelerator in the evolution of the 
galaxy. This would be a short-time-scale, one-time-only 
event as the galaxy approaches the group virial radius, 
where its star formation is quenched due to processes 
that are associated with the large gravitational potential 
of the massive group halo. A detailed study of depen- 
dence of star formation history of galaxies on these pa- 
rameters based the mass of the galaxy, and the local and 
global environments, will ultimately lead us to a better 
understanding of galaxy evolution. 

6. SUMMARY 

Using 905 galaxy groups identified by the pFoF algo- 
rithm in the RCSl photometric-redshift sample and cov- 
ering a group total halo mass range of ~ lO"'^^'^ M,- to 
~ IQ-'^^'^Mq, we study the group galaxy population over 
the rcdshift range 0.15 < z < 0.52. We examine the 
color-magnitude diagram of group galaxies as a function 
of redshift and group richness, and consider the effects of 
four parameters on the red galaxy fraction fj-ed and its 
evolution: the galaxy stellar mass M*, the total group 
stellar mass yi^^grp{as a proxy for the group halo mass), 
the group-centric radius rgrp, and local galaxy density 

We find the bright end (brighter than M^^) of the red- 
sequences in the CMD for stacked groups in redshift and 

group halo mass bins to be remarkably \miform, with 
their zero point, slope, and dispersion consistent with 
those found for clusters. Thus, the bright end of the 
red sequence is already in place, and likely formed at 
zTt2, even for those in groups approximately an order of 
magnitude less massive than clusters. 

Most of the evolutionary effects are seen in galaxies of 
lower stellar mass (strongest for M, slO-'^'^M©). We find 
that groups at lower redshifts possess larger /^ed than 
those at higher redshifts, exhibiting a group Butcher- 
Oemler effect. Examining the dependence of fred in more 
detail, we find: 



1. There is a strong dependence of fred on galaxy 
stellar mass. More massive galaxies have larger 
/red) and the group Butcher- Oemler effect is seen 
within our redshift range only for galaxies with 
M^^lOiiMg. 

2. The strength of the dependence of fred on the 
environmental parameters is also a strong hmc- 
tion of the galaxy stellar mass. Galaxies with 
M* felO^^M© are almost all red, independent of 
their local galaxy density, group-centric radius, 
and group halo mass. In the items that fol- 
low, the results apply primarily for galaxies with 

3. We find a dependence of fred on M^,^grp, in that 
galaxies in more massive groups have a larger fred- 
This is seen after the group galaxies are separated 
into Tgrp or E5 bins. This effect is strongest for the 
lowest mass galaxies. The group Butcher-Oemler 
effect appears to be stronger in lower-mass groups 
over this redshift range; i.e., the change of fred from 
the z ~ 0.5 to the 2; ~ 0.2 bin is larger for lower- 
mass groups. This difference is probably at least in 
part due to low-mass groups starting from a lower 
fred at 2; ~ 0.5. 

4. There is a dependence of fred on group-centric 
radius, which is stronger for lower-mass galaxies. 
This dependence still exists after controlling for lo- 
cal galaxy density, but is somewhat reduced. The 
fred-Tgrp trend is weak for the low-mass groups, 
especially in the two higher-redshift bins. While 
this difference may be real, part of it may be con- 
tributed by the larger uncertainty in determining 
the group center and/or larger contamination for 
the lower-mass groups. 

5. While group-centric radius and local galaxy den- 
sity are correlated, nevertheless, at a fixed Tgrp, 
there is still a significant dependence of fred on 
E5. The typical change in fred over a factor of 
about 13 in S5 is similar to that found in the lower- 
redshift SDSS "field" -galaxy sample, i.e., without 
controlling whether the galaxies are in clusters or 
groups. This change is also similar in magnitude 
to the change in fred seen in galaxies from outside 
the virial radius to the group core. The dependence 
on S5 is also larger for lower-mass galaxies. This 
indicates that group environment has a residual ef- 
fect over that of local galaxy environment (or vice 
versa), and the two must be considered at the same 
time. 

A general picture for galaxy evolution, and in partic- 
ular for galaxies in groups and clusters, emerges from 
these correlations and other work in the literature at dif- 
ferent redshifts. Within the four parameters that we have 
examined, in a 'nature versus nurture' scenario, we can 
consider the galaxy stellar mass as the predominant de- 
terminant of the evolutionary history of a galaxy. This 
would be the effect that produces the commonly accepted 
'down-sizing' of galaxy evolution. The environmental in- 
fluence on this secular trend, which accelerates the galaxy 
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to its final red state, can be considered as the 'nurtur- 
ing' part of the galaxy's history. There is sonic room left 
as to the debate of whether local galaxy density is part 
of the intrinsic character of the galaxy, in that galaxies 
born in a high density region may likely be preferentially 
more massive. However, local galaxy density, as an influ- 
ence on galaxy evolution may occur with different time 
scales and stochastically; e.g., galaxy harassment may 
be continuous, while galaxy mergers are a significantly 
shorter time scale event that occur more rarely. On the 
other hand, the eflfect of a galaxy falling into and being 
subsumed by a more massive halo is most likely a one- 
time-only effect that quenches the star formation in the 
galaxy over a relatively short time scale. These environ- 
mental events all have lasting effects on a given galaxy 
and the build-up of galaxy groups and clusters. To truly 
understand galaxy evolution, we need to have samples 
that cover a wide redshift range and are sufficiently large 
for us to clearly separate the effects of these parameters. 
This will then allow us to test the predictions, or inform 
the construction, of a variety of models and simulations. 



The recently completed RCS-2 (Gilbank et al. 2011) with 
a data set that will provide a similar sample of galaxy 
groups ~ 25 times larger than the current sample will 
allow us to make great strides toward this goal. 
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